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We study the role of the Hall current and electron inertia in collisionless magnetic reconnection within the 
framework of full two-fluid MHD. At spatial scales smaller than the electron inertial length, a topological 
change of magnetic field lines exclusively due to electron inertia becomes possible. Assuming stationary 
conditions, we derive a theoretical scaling for the reconnection rate, which is simply proportional to the 
Hall parameter. Using a pseudo-spectral code with no dissipative effects, our numerical results confirm this 
theoretical scaling. In particular, for a sequence of different Hall parameter values, our numerical results show 
that the width of the current sheet is independent of the Hall parameter while its thickness is of the order of 
the electron inertial length, thus confirming that the stationary reconnection rate is proportional to the Hall 
parameter. 


I. INTRODUCTION 

Magnetic reconnection is a physical process which 
converts magnetic free energy into kinetic energy and 
heat. This important mechanism of energy conversion 
is present in several space environments, such as so¬ 
lar flares and planetary magnetospheres^^— . The first 
model of magnetic reconnection was developed within 
the framework of one-fluid resistive magnetohydrody¬ 
namics (MHD), the so-called Sweet-Parker model^i^. In 
the Sweet-Parker regime, magnetic resistivity breaks the 
frozen-in condition at sufficiently small scales, thus al¬ 
lowing magnetic reconnection to occur. In particular, 
Parker— showed that the reconnection rate (i.e. the rate 
of change of magnetic flux due to reconnection) scales as 
the square root of the magnetic resistivity, which leads to 
exceedingly low reconnection rates for most space physics 
environments^. Years later, Petschek— reported a possi¬ 
ble way out to the slow-rate problem giving rise to the 
concept of fast magnetic reconnection, i.e. reconnection 
rates virtually independent of magnetic resistivity. In 
contrast to the Sweet-Parker scaling, the Petschek so¬ 
lution only showed a mild (logarithmic) dependence on 
magnetic resistivity, therefore being considered as fast 
reconnection. However, numerical results showed that 
the classical Petschek configuration cannot be attained 
in simulations with a spatially homogenous resistivity^. 

For collisionless magnetic reconnection, i.e. where the 
mean free path of the plasma particles is much larger 
than any length scale in the system (and therefore the 
resistivity is negligible), the resistive MHD model is 
no longer appropriate. From a more general perspec¬ 
tive, effects other than magnetic resistivity can break 
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the frozen-in condition, such as electron inertia or non- 
gyrotropic contributions to pressure. Discussions about 
the relative importance of electron inertia and non- 
gyrotropic (off-diagonal) pressure tensor terms can be 
found elsewhereiSr— . In the present paper we focus on 
the role of electron inertia. We also assume incompress¬ 
ibility and therefore non-gyrotropic pressure effects are 
not importanti More specifically, we focus on the phys¬ 
ical consequences of including the Hall effect and electron 
inertia (with isotropic pressures) into a fluidistic descrip¬ 
tion. 

At spatial scales larger than the ion inertial l ength Xj = 
c/ujpi (where c is the speed of light, Wpi = 4rKe^nlmi 
is the plasma proton frequency, e is the electron charge 
and n is the plasma density), the MHD description is ad¬ 
equate to describe global phenomena in most of the as- 
trophysical plasmas. However, at scales below Ai, where 
the ions become unmagnetized, the Hall-MHD (HMHD) 
description becomes valid. At spatial scales of the or¬ 
der of the electron inertial length Ae = c/wpe (where 
ujpe = sj 47re^n/me is the plasma electron frequency) or 
smaller, the terms of electron inertia become dominant, 
and electrons are no longer frozen to the magnetic field 
linesi. At this level of description, a topological change 
of the magnetic field lines exclusively due to electron in¬ 
ertia becomes possible. Andres et al.— presented a study 
of collisionless magnetic reconnection within the frame¬ 
work of Electron Inertia Hall-MHD (EIHMHD), i.e. a 
two-fluid theoretical framework that extends HMHD and 
includes the inertia of electrons. Using a pseudo-spectral 
code with no dissipative effects, the authors numerically 
confirmed that the change in the topology of the mag¬ 
netic field lines is exclusively due to the presence of elec¬ 
tron inertia. Moreover, they showed that the computed 
reconnection rates were independent of the mass ratio 
me/rui and remain a fair fraction of the Alfven velocity. 





2 


which therefore qualifies as fast reconnection. It is worth 
mentioning that the level of description of EIHMHD 
should not be confused with the so called electron MHD 
(EMHD) approximation. Instead, the EIHMHD model 
retains the whole dynamics of both the electron and ion 
flows throughout all the relevant spatial scales. It asymp¬ 
totically becomes MHD at the largest scales, HMHD at 
intermediate scales and EMHD at the smallest scales. 
Under the EMHD approximation, the ions are assumed 
to be static (because of their much larger mass) and the 
electrons are the ones to carry the electric current^. 

Geospace Environment Modeling (GEM) Reconnec¬ 
tion Ghallengeii was a project designed to study col¬ 
lisionless magnetic reconnection assuming different the¬ 
oretical approaches such as fully electromagnetic parti¬ 
cle in celU^^— , resistive MHD, HMHD^Si^— and hybrid 
codes^Sii^. The authors find that the reconnection rate 
is insensitive to the mechanism that breaks the frozen-in 
condition and its particular value is approximately ~0.1 
(in dimensionless form). In particular, Shay et al.— 
claimed that this values of the reconnection rate is an 
universal constant as the system become very large. 

However, several studies have demonstrated that 
the reconnection rate might still depend on the value 
of the Hall parameteri^S^— , the level of turbulent 
fluctuations^STi^ and the boundary conditions of the 
proble m^^i^^ . This idea that MHD turbulence may play 
an important role in a magnetic reconnection setup was 
first proposed by Matthaeus and Lamkin— . Smith 
et ai— examined the influence of the Hall effect and 
level of MHD turbulence on the reconnection rate in 2.5D 
compressible Hall MHD. Their results indicate that the 
reconnection rate is enhanced both by increasing the Hall 
parameter and by the turbulence amplitude. 

Following an approach of single-particle dynamics, 
Cowley— found an expression for the reconnection rate, 
which strongly depends on the ion inertia length. Wang 
et al.— reported an analytical treatment of quasi¬ 
stationary collisionless magnetic reconnection including 
the Hall effect, scalar electron pressure gradient and elec¬ 
tron inertia terms. The authors find that the reconnec¬ 
tion rate depends on the ion inertial length, the bound¬ 
ary/initial conditions and the expression for the external 
driving force. More recently and within the context of in¬ 
compressible HMHD, Simakov and Chacon— presented 
a quantitative analysis of reconnection valid for arbitrary 
values of the Hall parameter—. 

Our main goal in this paper is to study the magnetic 
reconnection rate, using a full two-fluid model for a com¬ 
pletely ionized hydrogen plasma, retaining the Hall cur¬ 
rent and electron inertia. Within this framework, we 
calculate a scaling for the quasi-stationary reconnection 
rate. Our results show that the reconnection rate has 
a linear dependence on the Hall parameter. In section 
m we briefly describe the ideal EIHMHD set of equa¬ 
tions. In section imi we present our theoretical scaling 
for the reconnection rate. In section IIV Al we show the 
set of equations that describes the dynamical evolution 


of the problem in a 2.5D setup and the corresponding 
initial conditions. In section HVBl considering a pseudo- 
spectral method to accurately run ideal simulations, we 
present our main numerical results. Finally, in section m 
we compare and discuss our results with those reported 
in the literature and summarize our main conclusions. 


II. ELECTRON INERTIA HALL-MHD MODEL 


The detailed derivation of the EIHMHD model have 
been presented elsewher o^^i^^ . In this section we sum¬ 
marize its key points. The equations of motion for an 
incompressible plasma made of ions and electrons with 
mass charge ±e, density rii = = n (because of 

quasi-neutrality), pressure pi^e and velocity respec¬ 
tively, can be written as 

rmn—^ = en(E -b x B) — Vpi (1) 

at c 

= —en(E -b x B) — Vpe (2) 
at c 


3 = 


Q 

— V X B = en(ui 
47r 


Ue) 


where 



clug 


dt 


+ (Ug.i • V)Ue,i 


( 3 ) 

( 4 ) 


is the total derivative. Here, B and E are the magnetic 
and electric fields, J is the electric current density and c is 
the speed of light. This set of equations can be written in 
a dimensionless form in terms of a typical length scale Lq , 
the constant particle density n, an intensity Bq for the 
magnetic field, a typical velocity va = B^j(3'KnMY^'^ 
(the Alfven velocity, where M = mi+rrie) and the electric 
field in units of Eq = vaBq/c, 

(1-m)^ = ^(E + WxB)-Vp, (5) 

M^ = -^(E + UgXB)-Vpg (6) 

J=^(W-Ue) (7) 

where we have introduced the dimensionless parameters 
/i = me/M and A = c/uJpmLq is the dimensionless Hall 
parameter, and ojpM = has the form of 

a plasma frequency for a particle of mass M. The di¬ 
mensionless ion and electron inertial lengths can be de¬ 
fined in terms of their corresponding plasma frequencies 
a}pi,e = simply as = c/wpi^eLe). Note 

that in the limit of electron inertia equal to zero, we ob¬ 
tain ujpM = Wpi, and therefore A = Ai = cfuipiLQ reduces 
to the usual Hall parameter. However, throughout this 
paper we are going to retain the effect of electron iner¬ 
tia through the parameter /i ^ 0. The expressions for 
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the dimensionless ion and electron inertial scales (Ai^e) 
in terms of the two dimensionless parameters /i and A 
are simply Ai = (1 — and Ae = /i^^^A. 

For a hydrodynamic description of this two-fluid 
plasma, we replace the velocity field for each species (i.e. 
Ui_e) in terms of two new vector fields. Namely, the hy¬ 
drodynamic velocity u given by 

U = (1 - n)Ui + flUe (8) 

and the electric current density J given by ©• From 
equations ©-®, we can readily obtain the velocity of 
each species as 

Ui = u + ^AJ (9) 

Ue = u — (1 — ^)AJ (10) 

The hydrodynamic equation of motion is the sum of the 
corresponding equations of motion ([5]) and m for each 
species 

^ = (11) 

where p = Pi +Pe is the total pressure. Even though most 
of the terms in equation (HU can easily be identified as a 
sum of the corresponding terms in equations ®-®, the 
sum of the convective derivatives in these equations are 
nonlinear terms that give rise to a new nonlinear term in 
equation (HU) which is proportional to p. Note also that 
in the limit of negligible electron inertia (i.e., for /i —>■ 0), 
equation (fTTll reduces to the equation of motion for the 
traditional one-fluid MHD. This is the case for the Hall- 
MHD description as well, which is also a two-fluid the¬ 
oretical description, but considering massless electrons 

{p = 0 ). 

On the other hand, the equation of motion for electrons 
(equation db])), using E = —dtA — S/cj) and (ug • V)ue = 
uje X Ue + V{u^/2) (with We = V X Ue being the electron 
vorticity) can be written as 

d 

— (A - p\Ue) = Ue X (B - pXUe) + 

+ {Xpe + pX-^ — 4^) ( 12 ) 

We define, 

B' = B — pXue = B — p(l — p)}?V^B — pXu) (13) 

where w = V x u is the hydrodynamic vorticity. Taking 
the curl of equation (1121) we obtain a dynamical equation 
for the magnetic field 

a* B'= V X [u-(l-/r)AJ] X B' (14) 

Equations HU and (fT^ are the EIHMHD equations. It 
is interesting to note that the presence of the electron 
mass introduces higher order derivative terms. This cer¬ 
tainly has an impact at large wavenumbers, affecting the 
distribution of energy at very small scales. Note that in 
the limit of negligible electron inertia (i.e., for p ^ 0), 
equations HU and (fTTl) reduce to the standard equation 
of motion and induction equation of 



Figure 1. Schematic 2.5D reconnection region. 


III. THEORETICAL SCALING OF THE 
MAGNETIC RECONNECTION RATE 


In the context of collisionless magnetic reconnection, 
the reconnection region develops a multi-scale structure 
in which the ion and electron inertial lengths Xi^e play 
a rolei^. As we discussed in the Introduction, ions can 
be considered approximately static and electrons are the 
ones to carry most of the electric current. Also, at these 
scales the terms of electron inertia become dominant, 
and the electrons can no longer be frozen-in to the mag¬ 
netic field linesi. Therefore, at this level of description, a 
change in the topology of the magnetic field lines which 
is exclusively due to electron inertia, becomes possible. 

Within scales near the X-point, where ju^j << |J/en|, 
we obtain a scaling for the reconnection rate as a func¬ 
tion of A and p which are the main parameters of the 
problem. We consider a rectangular reconnection region 
with a width 2S and a length 2 A (see Figure [U- By defi¬ 
nition, the reconnection rate in a 2D configuration is the 
out-of-plane component of the electric field (i.e., E^) at 
the X-point. The electric field can be obtained from the 
ideal equation of motion for the electrons dU) as 



, e 

+ Ue X (We H- B) 

at rrieC 


+V 



Pe 

meTL 


(15) 


Under the assumption of quasi-stationarity (i.e., dt 0) 
for a 2.5D setup (i.e., dz ^ 0), the out-of-plane compo¬ 
nent of the electric field (the z direction) reduces to 


Ez = -^ Z • Ue X We = - 5 -^ Z • J X (V X J) (16) 

e 


where we have assumed Ug ~ —Ijen. 

In view of the sketch shown in Figurejl] close to the X- 
point is dx ~ A“^, dy ^ and Jz = cBinjAnb, where 
Bin is the magnetic field at the edge of the reconnection 
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region in the inflow direction. Therefore, 

/ C \ ^ B^Bin 
^ e \4TTne) 


(17) 


To estimate the out-of-plane component of the magnetic 
field we consider the z component of the curl of 

equation (fT51) (under quasi-stationary conditions), i.e. 


J X ( -^B - — V X J 

CTTle en 


z • V X 

which, in 2.5D setup leads to 


= 0 


w. 


pe 


(18) 


Bu • Vi J, = - Ji • Vi(V2B,) (19) 


and therefore 


V • u = 0, we can write the magnetic and velocity fields 
as 


B = V X [z a{x,y,t)] + z b{x,y,t) (24) 

u^V x[z ip{x,y,t)]+zu{x,y,t) (25) 

where a{x,y,t) and Lp(x,y,t) are the scalar potential for 
the magnetic and velocity fields respectively and b{x, y, t) 
and u{x,y,t) are simply the corresponding out-of-plane 
components. In terms of these scalar potentials, equa¬ 
tions CD and (fT^ take the form 

dt uj= [ip, uj] - [a, j] - (1 - m)M^ [b, V^fe] (26) 

dtu = [p, m] - [a, 5] - (1 - y)n\^ [j, b] (27) 

dt a' = [ip — {1 — /i)A6, a'\ (28) 

dt b' = [p- {I- fi)Xb, b'] -I- [m - (1 - fj.)Xj, o'] (29) 

where 


B, = ( 20 ) 

c 

The z-component of the electric held at the X-point is 
then 


oj = — 

(30) 

j = -V^a 

(31) 

a' = a + {1 — p)p}?j — pXu 

(32) 

b' = b — {1 — p)pX^V^b — pXu) 

(33) 


= c 

Airen ASujpe 


( 21 ) 


The dimensionless reconnection rate, i.e. r = cEz(Bqva^ 
becomes 


r = 


c c 



( 22 ) 


As it was discussed in the Introduction, we expect Bin 
and A not to depend on Their particular values are 
only determined by the boundary and initial conditions. 
Nevertheless, in the next section we evaluate the poten¬ 
tial dependence of Bin , ^ and A with the Hall parameter 
in our numerical results. 

Assuming that the thickness of the current sheet is 
essentially the electron inertial length, i.e. 6 ~ cjupe and 
also that the typical magnetic held intensity is Ho = Bin 
and the typical length scale is Lq = A we obtain 


r = X (23) 

Note that if 5 ^ cjupe, according to (l20)) we also ob¬ 
tain that (in the regime of quasi-stationary reconnection) 
Bz ~ Bin- Note also that the reconnection rate is inde¬ 
pendent of the mass ratio as shown in Andres et al.—. 


IV. NUMERICAL RESULTS 

A. 2.5D Setup and initial conditions 

In a 2.5D setup, the vector helds depend on two co¬ 
ordinates, say X and y, although they have their three 
components. Considering the incompressible case, i.e. 


and the nonlinear terms are the standard Poisson brack¬ 
ets, i.e. [p, g] = dxpdyq — dypdxQ- The set of equations 
(l26l) - (l2^ describe the dynamical evolution of the mag¬ 
netic and velocity helds in 2.5D. When p = 0 (massless 
electrons) this set of equations reduces to the incompress¬ 
ible 2.5D HMHD equations^. 

In the present paper, we performed 2.5D 
EIHMHD simulations using a pseudo-spectral code, 
which yields exponentially fast numerical convergence 
and negligible numerical dissipation. The accuracy 
of the numerical scheme can be verihed in part by 
looking at the behavior of the ideal invariants of the 
EIHMHD equations in time. The simulations reported 
here correspond to zero viscosity and resistivity, and the 
total energ y^ is conserved by the numerical scheme 
with an error AE/E of less than 10“®. The ion and 
electron helicities were initially zero, and throughout 
their evolution differ from zero in less than 10“^®. 
Therefore, hereafter we assume that our code conserves 
energy. The reconnection processes that are observed 
to occur, must then be the exclusive result of electron 
inertia. 

Our initial condition to simulate a thin current sheet 
is given by (assuming periodic boundary conditions in a 
27r X 27r box) 


B{x,y,t = 0) = Bo 


tanh 



— tanh 



(34) 


where, in normalized units, we have Bq = 1 and I = 0.02. 
To drive reconnection, a monochromatic perturbation 
(5B = V X [z 6a{x, ?/)] with 6a{x, y) = Oq cos^k^x), fca, = 1 
and an amplitude of oq = 0.02Bo is added to the initial 
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Figure 2. The image (in grayscale) shows the spatial distri¬ 
bution of current density j{x,y) at t = 1.0 for A = 0.1 and 
rrie/mi = 0.015. Contour levels of a{x,y) are superimposed 
(black lines). 

condition (1341) . It is worth mentioning that in our simula¬ 
tions we do not use any external driving force, and there¬ 
fore the reconnection process can be regarded as self- 
driven. We perform numerical simulations with a spatial 
resolution of 2048^ grid points. For all the runs we use a 
value of electron to proton mass ratio me/wi = 0.015 and 
different values of the Hall parameter A. Figure [5] shows 
the set up of magnetic reconnection for A = 0.1. Contour 
levels of magnetic flux a{x, y) are in black lines, super¬ 
imposed to the electric current density component along 
the z direction, j{x,y), at time t = 0.6 (in grayscale). 
We only show half a box of integration for each case, of 
size 27r x tt. 


B. Quasi-stationary magnetic reconnection 

Within the framework of EIHMHD, we study the col¬ 
lisionless magnetic reconnection problem varying the di¬ 
mensionless Hall parameter A. Using the initial con¬ 
ditions described in subsection IIV A1 we performed ten 
ideal runs with a spatial resolution of 2048^ grid points 
for different values of the Hall parameter. Our runs 
span the range A = 0.07 to A = 0.16, with a step of 
0.01. The values of A are sufficiently small, to minimize 
the potential influence of boundary conditions. In all 
these runs the electron to ion mass ratio corresponds to 
me/rrii = 0.015. 

To measure the efficiency of the magnetic reconnec¬ 
tion process, the dimensionless reconnection rate r(t) is 
dehned, which is the rate at which magnetic flux flows 
into the X-point. Using equation (0^ it is straightfor¬ 
ward to show that the total reconnected flux <I>(t) is 
4)(t) = amax—amirr^^- Therefore, the reconnection rate 
r{t) is the variation of the magnetic flux per unit time, 
i.e. r{t) = d^{t)/dt. Figure [3] shows the reconnected flux 
(upper panel) and reconnection rate (lower panel) as a 
function of time, for the ten values of the Hall parame¬ 
ter. In contrast to previous claim s Figure |3] shows 
that the reconnection rate strongly depends on A and is 
not a universal constant. As it can be seen, the recon¬ 
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Figure 3. Reconnected flux 4> (upper panel) and reconnec¬ 
tion rate r (lower panel) as a function of time for A = 
0.07, •• • ,0.16 (from bottom to top). For all runs the elec¬ 
tron to ion mass ratio is m^/mi = 0.015. 


nected flux monotonically increases with A, which ulti¬ 
mately leads to an increment of the maximum magnitude 
of reconnection rate. We also note that as we increase A, 
the maximum reconnection rate occurs at earlier times. 
Similar behavior has been reported in the literature when 
the Hall effect is included in Ohm’s la w^^T^T^ . 

From equation (|2^ we see the importance of studying 
whether the thickness and length of the reconnection re¬ 
gion ((5 and A, respectively) and the magnetic field at 
the edge of this region {Bin) change as a function of 
the Hall parameter. Since our scaling was performed as¬ 
suming quasi-stationary conditions, we have to take this 
constraint into account. The width of the reconnection 
region 5 is dehned in terms of the current density pro¬ 
fit® j{v) across the layer—. The value of 5 is obtained 
from a best ht of the numerical prohle to a sedc?{y/5) 
function, which is consistent with the initial prohle give 
by equation (I34|) . To determine Bin we simply adopt 
Bin = Bx{x = 'Kl2,y = 7r/4 — <5), since our neutral 
point is located at x = 7r/2,2/ = 7r/4. We assume that 
the system evolves in a quasi-stationary fashion during 
a time interval such that 5 and Bin show approximately 
no temporal variations. The length of the reconnection 
region A was obtained from the outhow velocity prohle 
Ux{x, y = 7r/4) applying the incompressible condition for 
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Figure 4. Quasi-stationary values of 5 (gray circles) as a func¬ 
tion of Ae- We plot the electron inertial length Ae in gray- 
dashed line for reference. 


Figure 6. Quasi-stationary reconnection rate r (gray circles) 
as a function of the Hall parameter A. The best linear-fit for 
log A—log r is shown in gray-dashed line. Inset: Ratio between 
quasi-stationary reconnection rates and the Hall parameter 
(gray squares) as a function of the Hall parameter. 
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Figure 5. Quasi-stationary values of Bin (upper panel) and 
A (lower panel) as a function of A. The gray-dashed line 
indicates the mean values of Bin and A. 


the plasma, i.e. 


= I + A,y = J) = ^ 4“)(x =^,y = ^-S). 


(35) 

Figure S] shows the quasi-stationary values of S (gray 
circles) as a function of Ag. In addition, we plot Ag in 
gray-dashed line. As expected, the width of the reconnec¬ 


tion region is of the order of the electron inertial length. 
In particular, from a best linear-fit for log 6 — log Ag we 
obtain S = (1.3 ±0.3) Ag °®^° *’^. Therefore, we conclude 
that 5 ^ Ag. 

Figure [S] shows Bin (upper panel) and A (lower panel) 
as a function of A (gray squares) for the ten values of the 
Hall parameter. Figure [5] indicates that Bin and A show 
approximately no dependence with the Hall parameter. 
This result is compatible with previous results reported 
in the literatur o^^i^^ . 

The results displayed in Figure 0] ((5 ^ Ag) and Figure 
El (Bin ^ const and A ~ const) lend support to the as¬ 
sumptions made in equation (l22ll to obtain equation (1231) , 
i.e. that the reconnection rate is simply proportional to 
the Hall parameter. Figure [6] shows the quasi-stationary 
reconnection rates (gray circles), i.e. the mean recon¬ 
nection rate for the time interval determined in Section 
II V B1 as a function of the Hall parameter A. In addition, 
we plot the curves corresponding to the best linear-fit for 
log A—logr (dashed line). The inset in FigurelHlshows r/A 
(gray squares) as a function of A. From the best linear-fit 
for logA - logr we obtain r = (0.11 ± 0.07) 

Therefore, we conclude that the reconnection rate r is 
compatible with a linear relation with the Hall parame¬ 
ter A, as it was predicted by our analytical relation (1231) . 

Finally, we also compare the quasi-stationary recon¬ 
nection rate for a fixed value of the Hall parameter (A = 
0.1), and two different electron to proton mass ratios. In 
particular, we compared the results for Wg/mi = 0.015 
and rUg/TOi = 0.15. In the quasi-stationary regime, we 
find approximately the same reconnection rate. This re- 
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suit is compatible with our theoretical result, which pre¬ 
dicts that fast reconnection rate is insensitive to the elec¬ 
tron to proton mass ratio even though it needs to be 
nonzero for reconnection to take place (Birn et al., 2001; 
Zenitani et ai, 2011; see also Andres et ai, 2014a). 


V. DISCUSSION AND CONCLUSIONS 

Within the framework of two-fluid MHD and assuming 
stationary conditions, we obtain a theoretical scaling for 
the reconnection rate. Our numerical results confirm our 
assumptions that the thickness of the current sheet is 
essentially the electron inertial length, i.e. d ~ Ag, and 
that Bin and A do not depend on the Hall parameter—. 
More importantly, our numerical results also confirm the 
predicted linear dependence of the reconnection rate r 
with the Hall parameter A (i.e. r oc A). 

Within the context of incompressible HMHD, Simakov 
and Chacon— presented a quantitative analysis of re¬ 
connection valid for the resistive, HMHD and EMHD 
regimes. Their study concentrated on the reconnection 
region, without considering any particular external driv¬ 
ing force. In the resistive MHD limit, the authors recover 
the standard resistive result^. In the limit of EMHD, the 
authors find that the reconnection rate does not explic¬ 
itly depend on the dissipation coefficients and features a 
strong dependence on the Hall parameter. In particular, 
they confirm an earlier result and hnd that r = v^A/Z^— , 
which is consistent with our scaling. 

Malyshkin— also calculated the rate of quasi¬ 
stationary, 2.5D magnetic reconnection within the frame¬ 
work of incompressible HMHD. The author find that the 
dimensionless reconnection rate is independent of the 
electrical resistivity and equal to A/L, where L is the 
scale length of the external magnetic field in the upstream 
region outside the electron layer. This result is also com¬ 
patible with our theoretical results^. 

In a different direction, Wang et al.— reported a sim¬ 
ilar linear dependence with A and noted that Bm is de¬ 
termined by the functional form of the boundary con¬ 
ditions, while A depends on a external time-dependent 
driving force. For a particular model of external driving, 
Wang et ai— calculated the scaling of the reconnection 
rate within the framework of resistive HMHD. The au¬ 
thors found a A^/^ dependence for the reconnection rate. 
This particular scaling is not comparable with our re¬ 
sults, since in our simulations we do not consider any 
external driving force. 

As discussed in the Introduction, MHD turbulence 
may play an important role in magnetic reconnection^^. 
Smith et ai— examined the influence of the Hall effect 
and level of MHD turbulence on the reconnection rate 
in 2.5D compressible Hall MHD. Their results indicate 
that the reconnection rate is enhanced both by increasing 
the Hall parameter and by the turbulence amplitude. In 
agreement with these studies, our numerical results show 
a clear enhancement as we increase the Hall parameter. 


Smith et ai— also suggested a power-law scaling of the 
reconnection rate as a function of the Hall parameter as 
r oc 5B\^/^, where SB is the level of initial turbulence 
in the system. However, in our study we do not consider 
any initial turbulence level, since we focus on the con¬ 
sequences of adding the Hall effect and electron inertia 
terms in a laminar background. Also, in their simula¬ 
tions, Smith et al— added a small amount of magnetic 
resistivity, in order to break the frozen-in condition and 
start the reconnection process, which is different from our 
ideal EIHMHD description. 

In summary, we obtained a theoretical linear scal¬ 
ing for the reconnection rate as a function of the Hall 
parameter, which is confirmed by our numerical re¬ 
sults and is also compatible with previous results in the 

literatur o^^i^^ 
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